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ABSTRACT 



The computerized measurement of thermoacoustically generated 
temperature gradients in short, thin plates is reported. The 
computerized data acquisition system is delineated. The 
temperature difference developed across a stack of short plates 
was measured as a function of the longitudinal position of the 
plates in a resonant tube for acoustic pressure amplitudes of 0.5 
to 6.6 kPa, and static (or mean) pressures from 100 to 440 kPa, 
in argon and helium for the first through the third harmonic 
frequencies of the tube. Measured data were compared with 
predictions based on work done by Wheatley and others [J. 
Wheatley, et al.. Journal of the Acoustical Society of America, 

V. 74, pp. 153-170, 1983] and results reported by Muzzerall 
(Master's Thesis in Engineering Acoustics, Naval Postgraduate 
School, Monterey, CA, September 1987). For low acoustic and 
static pressures, there is good agreement between measured data 
and theory. As the acoustic pressure amplitudes increase there 
is a general degradation of agreement up to the point at which it 
appears saturation of the thermoacoustic effect occurs. 
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I. 



INTRODUCTION 



The investigation of thermoacoustic heat engines was begun by 
Wheatley and associates in the early 1980's [References 1 & 2] . 
They explored the phenomena upon which these engines are based by 
introducing a thermoacoustic couple (TAG) into an acoustic 
standing wave. A TAG is a stack of short, thin plates. At least 
one of those plates is instrumented with a series of 
thermocouples to measure temperature difference developed across 
it as a result of the heat transport. In 1987, Muzzerall [Ref. 

3] made measurements of the temperature difference developed 
across a TAG as a function of the position of the TAG in the 
acoustic field and compared the results to the theoretical value 
given by Wheatley, et al. There was good agreement between the 
measured data and the theory at low acoustic pressure amplitude 
however, at high pressure amplitude there was agreement only in 
the regions near the pressure antinodes. Elsewhere, 
unaccountable severe skewing and an overall reduction in the 
temperature difference was observed. 

Muzzerall published extensive data only at two acoustic 
pressure amplitudes separated by 12 dB. Also, his spatial 
resolution was only 1 cm. The purpose of the experiment reported 
in this thesis is to extend this investigation to see where the 
theory and measured data start to diverge and to further explore 
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what happens at high acoustic and mean pressures. In order to 
accomplish these goals, better spatial resolution and the ability 
to investigate a broader range of parameters is required. The 
parameters to be controlled are acoustic pressure amplitude, 
ambient (or mean) pressure, frequency, and gas type. In order to 
handle the quantity of data recorded and to ensure repeatability 
of measurements, both the data acquisition and the linear 
positioning of the TAG are computer controlled in this 
investigation. Following a derivation of the theoretical basis 
for the generation of temperature gradients and a summary of 
Muzzerall's results, the computer controlled experimental 
apparatus is delineated. Methodology for data acquisition is 
discussed and a graphical comparison of measured and theoretical 
data is made. Areas for further investigation are included. 
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II . THEORY 



Wheatley, et al. [Ref. 2] used the term "natural engine" to 
describe a process where a phase shift is introduced between 
temperature and velocity by a natural, irreversible process - 
here thermal conduction. This process leads to, what has been 
termed, thermoacoustic heat transport. The TAG acts as such an 
engine establishing a thermal lag which drives an entropy flow. 
A brief description of this phenomena is described below. A 
detailed development of theoretical thermoacoustics applied to 
engines is given by Swift [Ref. 4] . 

Following the description of the phenomenon, a summary of 
Muzzerall's findings is presented. 



A. BASIC THEORY 



The cyclic heat pumping action of the TAG as a "natural 
engine" may best be described from a Lagrangian point of view. 
Referring to Figure 1, a poor thermally conducting plate is 
located near the end of a rigid resonant tube. A sound source is 
located at the left end of the tube (not shown). This source 
generates an acoustic standing wave in the tube. Consider a 
parcel of gas located near (within a thermal penetration depth) 
the plate. Both the plate and the parcel of gas are initially at 
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Rigid Closed Tube 




Figure 1 

The Thermoacoustic Effect 
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temperature T. During the compression phase of the acoustic 
cycle the parcel is displaced toward the pressure antinode and 
undergoes an adiabatic compression increasing it's temperature to 
T++. The temperature difference drives a heat flow (Q) from the 
gas to the plate reducing the parcel's temperature to T+. On the 
Expansion phase of the acoustic cycle, the gas undergoes an 
adiabatic expansion, but is now at the lower temperature T-. 

Here the temperature difference results in a heat flow in the 
opposite direction, returning the gas parcel to its original 
temperature T. The result is that after one cycle, the parcel 
has transported an amount of heat Q over a short distance. The 
deposited heat is picked up by an adjacent parcel and moved still 
further during the next cycle. The heat is shuttled down the 
plate, toward the closed end of the tube, until it reaches the 
end of the plate. There it collects, until conducted back 
through the plate or the surrounding gas. 

The system reaches steady state when the hydrodynamic heat 
(entropy) flow in the gas is balanced by the return diffusive 
heat (entropy) flow in the plate and in the surrounding gas. The 
entropy flow increases when there is no average temperature 
gradient in the plate. As the gradient increases, the heat flow 
begins to drop. There exists a critical gradient beyond which 
acoustics alone will no longer increase the temperature gradient 
in the plate [Ref. 1] . The gradients from zero through the 
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critical value demark the heat pump regime, which is addressed 
specifically in this thesis. 



The behavior of the temperature difference as a function of 
its position in the standing wave can be derived from Wheatley’s 
estimate for the time-average hydrodynamic heat flow [Ref. 1: p. 
156, Eqn. 15]: 



Q-7T 5i. 

where 

TT = overall perimeter of the plate 
= thermal penetration depth 
= mean temperature 
Q = isobaric expansion coefficient 
Pj^ = acoustic pressure magnitude 
Uj^ = acoustic velocity magnitude 
P = defined by Wheatley as the ratio of 
The terms I’m/? determine the adiabatic temperature change per 
pressure change and 7T5 ^ is the volume of gas-plate interaction. 
(P- 1) shows how the entropy flow changes as the temperature 
gradient changes and becomes zero at Vt = Vt^,j.^^. The next 
step in the development [Ref. 1: p. 158, Eqn. 30] is to note 
that : 
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At = <Q2>/(kA/l) 

where 

k = thermal conductivity 
A = sectional area of the plate 
1 = effective length of couple 

<□2 > = spatial average with subscript 2 denoting second order 
validity 

The relationship between acoustic pressure, velocity and the 
temperature difference is seen in Figure 2 for arbitrary units. 



Wheatley's expression for the temperature difference across 
the TAC [Ref. 2: p. 157, Eqn. 17] is: 



At = 



K 5k (1 sin(2KX) 




Pq = acoustic pressure (AC pressure) 
O' = Prandtl # 

= thermal penetration depth 
a = sound speed 
K = wave number 
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Figure 2 

Spatial Dependence of Pressure, Velocity and T 
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X 



distance from closed end of tube 



= mean density 

kp , kg = thermal conductivity for plate(p) and gas(g) 

1 = length of TAC 
Tjj, = TAC temperature 
CJ = 27T . frequency 

ratio of specific heats, Cp /C^ 

which is what was used to determine the theoretical values in 
this paper, 

B. MUZZERALL'S RESULTS 



Included for comparison are Muzzerall's results (Figures 3 
and 4) for a TAC made of stainless steel plates 2.26 cm long x 
2.16 cm wide, at four penetration depth spacing. These 
measurements were made in argon at a frequency of 705Hz and a 
static pressure of one atmosphere. The lines indicated as theory 
are the predictions of the expression given previously for ^T. 
The acoustic pressure amplitude in Figure 3 is 630 kPa. There is 
overall good agreement between theory and experiment. In Figure 
4 the acoustic pressure amplitude has been increased to 2500 kPa. 
The overall agreement between theory and experiment is severely 
diminished, although it is still quite good in regions near 
acoustic velocity nodes. The impetus behind the present thesis 
is to explore the regions between and beyond these two figures in 
order to better understand the reasons for this behavior. 
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Figure 3 

Muzzerall’s Results 150dB 
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Figure 4 

Muzzerall’s Results 162dB 
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III. APPARATUS 



The experimental apparatus can most easily be described by 
breaking it down into its functional component parts. The 
subsystems are: the tube and pressure vessel, the TAC; the 

linear positioning system; and the data acquisition (DACQ) 
system. A brief description of the function of these subsystems 
is given along with block diagrams for each. 

A. TUBE AND PRESSURE VESSEL 

1. Resonant Tube 

The tube is a 1.22 m long, 3.81 cm inside diameter 
copper tube. The length was chosen to obtain one half of an 
acoustic wavelength - therefore a full thermal wavelength - at 
the fundamental. The tube is soft soldered to brass flanges 
which are attached to the pressure vessel at one end and a 
blanking plate at the other. The purpose of the blanking plate 
is to provide a rigid termination. An Endevco model 8510B-5 
piezoresistive transducer is threaded through the blanking plate 
and used to measure the acoustic pressure. A 1/8" diameter 
stainless steel probe tube, which carries the TAC, also goes 
through the blanking plate. It is supported in the center of the 
plate by a rubber 0-ring compression fitting. In order to allow 



12 



for elevated mean pressures within the tube, the volume behind 
the Endevco piezoresistive diaphragm is connected to the resonant 
tube via a capillary leak, which allows pressure equalization, 
but blocks acoustic pressure oscillations. The tube is wrapped 
with insulation material to ensure thermal isolation. 

( Figure 5 ) . 

2. Pressure Vessel 

The ability to increase mean pressure required the 
acoustic driver (discussed below) to be enclosed in a pressure 
vessel. The vessel is made out of an aluminum cylinder joined to 
two aluminum end plates by threaded rods. The cylinder is 30.5 
cm long with a 25.4 cm inside diameter and 2.54 cm wall 
thickness. The end plates are 2.54 cm thick disks, 35.6 cm in 
diameter. The mean pressure inside the pressure vessel is sensed 
both by a dial pressure gauge and by an OMEGA Model PX304-150AV 
pressure transducer mounted in the back end plate. A 150 psi 
relief valve is mounted on the pressure vessel end plate. 

3 . Driver 

The acoustic driver is a JBL Model 2445J compression 
driver with input provided by a Hewlett Packard Model 3314A 
Function Generator amplified by a Techron Model 7520 Power Supply 
Amplifier. The driver is mounted on the inside of the pressure 
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Figure 5 

Data Acquisition System 
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vessel end plate attached to the resonant tube. Provisions are 
made to enable mounting a positive displacement piston driver for 
work at high ambient pressures, though it was not used in this 
experiment . 

4. Plumbing 

Connections were made to pressurize the tube and 
pressure vessel, with provisions to equalize the pressure on both 
sides of the acoustic driver. Connections were also made to 
attach a vacuum pump, with a needle valve for controlled 
backfilling, and also a system vent for outlet. 



B. TAC 

The TAC is a stack of 6 G-10 fiberglass plates, 1 cm long by 
2.54 cm wide (Figure 6). The center two plates are glued 
together with four type E ( chromel-constantan ) thermocouples in 
between. The thermocouples are connected in series to increase 
the output voltage and they allow the measurement of temperature 
difference across the TAC. The spacing between the center 
"sandwich" and the remaining guard plates is 1.54 mm. A 
thermocouple is glued to the back edge of one of the guard 
plates, which allowed the measurement of the absolute TAC 
temperature. The center plate sandwich is 0.595 mm thick and the 
guard plates are 0.36 mm thick. 
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C. DACQ 



Data acquisition is centered around the HP 3457A multimeter 
(Figure 5). To reduce the number of connection lines in the 
figure, circled numbers represent connections to multimeter input 
channels and the circled "f" the front panel input. The Kikusui 
Model 6100 Oscilloscope was used to monitor system performance. 

The system reference temperature was read from a thermistor 
which was placed inside an aluminum block to ensure thermal 
stability (Figure 7). Junctions for tube and TAG temperature 
thermocouple wires are attached to this block, insuring the 
chromel-copper and constantan-copper junctions are isothermal. 



D. LINEAR POSITIONING 



Motive force for TAG positioning is provided by a Gompumotor 
Model M83-135 Stepper Motor (25,000 steps/rev) (Figure 8) which 
is coupled to a Warner Electric 0.125 Lead Screw (8 turns/inch) 
on which rides a precision ball nut having zero lash. The ball 
nut is attached to an arm which slides on a support rod for 
stability. The other end of the arm is attached to the stainless 
steel probe tube which holds the TAG. A Gompumotor Model 2100 
Indexer serves as the interface between a Standard 286 (AT clone) 
computer and the motor . 
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Figure 7 

Temperature Reference 
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Figure 8 

Linear Positioning System 
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AT CLONE GPIB 



End of travel limit for both CW and CCW rotation and 



mechanical HOME (or reference) limit switches were tied directly 
to the indexer. All data runs start from the HOME position. To 
more precisely define HOME, a 1000:1 optical encoder is used in 
conjunction with the mechanical home limit switch. Once the 
mechanical switch activates, rotation is continued until the 
appropriate index on the encoder wheel is found. Positional 
repeatability was consistently within 1 motor step, or 
approximately 0.127 microns of longitudinal travel, assuming the 
mechanical switch was consistent enough to bring the optical 
encoder into the proper range. 

Interface between the personal computer (PC), the indexer, 
and the other HP equipment is accomplished through a National 
Instruments IEEE 488 General Purpose Interface Bus (GPIB) card in 
the PC. 
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IV. MEASUREMENT 



Development and testing of the experimental apparatus with 
its complex control and data collection program was a major time 
consuming undertaking. However the result is that once the 
system is manually pressurized, all subsequent operations are 
fully automated. 

Data acquisition is highly flexible. For example, the 
ability to adjust linear positioning step size ensured a data 
density sufficient to show the most detailed features of the 
temperature profile. Experimental results are presented in a 
serial manner with one feature changing per series. A summary 
figure leads to the discussion section. 



A. PROCEDURE 

Before data acquisition is started, the system is pumped down 
and backfilled with either helium or argon several times. 
Approximate resonance frequencies are determined next by driving 
the JBL with the noise source of the Signal Analyzer. The 
Endevco output is then fed back into the Signal Analyzer. The 
displayed signal peaks denote resonance frequencies. The driver 
input device is now changed to the function generator which is 
set to individual resonance frequencies to check for interference 
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between those resonances and higher harmonics of the drive signal 
(Figure 9). In general the harmonics were 18 or more dB below 
the fundamental . 

The gain on the voltage amplifier was set to a mid-range 
position and then left in that position for the duration of the 
data run. Determination of the magnitude of the amplification 
was unnecessary as acoustic pressure was measured directly. 

Automated data acquisition now begins. The program initially 
takes voltages from the ambient pressure transducer and directly 
converts them to pressure. The reference temperature is then 
read directly from the multimeter and converted to a effective 
reference voltage through a polynominal. The low order 
polynomial was numerically determined from NBS thermocouple 
tables over the limited temperature range of interest. Voltages 
from both the tube and TAC temperature thermocouples are read, 
added to the effective reference voltage, and converted to 
temperatures by an inverse polynominal. The drive voltage 
amplitude is set on the function generator and a search is made 
for resonance frequency by stepping on both sides of the 
previously determined natural frequency and reading the resultant 
acoustic pressure amplitude. The frequency producing the highest 
acoustic pressure amplitude is selected and set in the function 
generator. 
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RANGE; 5 dBV STATUS: PAUSED 

A: MAG RMS: 20 




X; 216.875 Hz Y: 5.51 dBV 



Figure 9 

Drive Signal Spectrum 
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Once the preliminary readings are finished and the function 



generator set, the data run commences by reading the TAG 
thermocouple voltage and converting it to a temperature as 
discussed above. Then the TAG thermopile voltage is read and 
converted to a temperature differences as follows: 



dV = TAG thermocouple voltage 
dT = TAG temperature 

The TAG is then moved to the next position, a wait time, obtained 
from Muzzerall [Ref. 3], activated and then readings are taken. 
Multiple data runs were made at constant static pressure and 
increasing acoustic pressure by resetting the indexer to the home 
position and incrementing the drive voltage on the function 
generator. 

B. RESULTS 

Experimental results are presented here in several series 
interspersed with theoretical predictions. The ordinate on all 
of the figures is the temperature differences (T„ - T^, ) in 
degrees Kelvin across the TAG. The abscissa is the product of 
the wave number and the distance of the center of the TAG from 
the rigid end of the resonance tube. All pressures, acoustic and 




Av dT ( type E ) 



dV 



J T AC T 




TAG thermopile V 
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ambient are in kilopascals. The first plot (Figure 10) is the 
only one to show individual data points - all other data is shown 
by connecting lines. The purpose of this figure is to show the 
density of data points, which is similar for all the subsequent 
figures . 

The first series (Figures 11, 12, and 13) demonstrate the 
typical degradation of experimental results against theoretical 
prediction at increasing acoustic pressures. The next series 
(Figures 14 through 17) show the effect of shifting from the 
first through the third overtone. Remaining at the third 
overtone, the next series (Figures 17 through 21) show increasing 
ambient pressure in one atmosphere steps. The last series 
(Figures 22 through 26 ) shifts to Helium at one atmosphere and 
shows again increasing frequencies. 



C. DISCUSSION 

The graphical presentation of the "raw" experimental results 
show no sharp demarcation between agreement of theoretical 
prediction and measured data, but demonstrate a gradual degration 
of agreement at increasing frequencies and pressures. Several 
interesting features can be seen in the series. 

The first series demonstrate results akin those of Muzzerall 
[Ref. 3]. Specifically shown is good agreement at velocity nodes 
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Figure 10 

Individual Data Points Plotted 
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ARGON 416 Hz 1.328 kPa 330 kPa 
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Figure 11 

Series la Increasing Acoustic Pressure 
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Figure 12 

Series lb Increasing Acoustic Pressure 
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ARGON 415 Hz 3.184 kPa 330 kPa 
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Figure 13 

Series Ic Increasing Acoustic Pressure 
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Figure 14 

Series 2a Increasing Frequency 



30 



ARGON EXPERIMENTAL 217 Hz 103 kPa 
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Figure 15 

Series 2b Increasing Frequency 
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Figure 16 

Series 2c Increasing Frequency 
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Figure 17 

Series 2d Increasing Frequency 
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Figure 18 

Series 3a Increasing Static Pressure 
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ARGON THEORETICAL 415 Hr 330 kPa 
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Figure 19 

Series 3b Increasing Static Pressure 
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ARGON EXPERIMENTAL 415 HZ 330 kPa 
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Figure 20 

Series 3c Increasing Static Pressure 
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Figure 21 

Series 3d Increasing Static Pressure 
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HELIUM THEORETICAL 694 Hz 107 kPa 
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Figure 22 

Series 4a Increasing Frequency 
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HELIUM EXPERIMENTAL 694 Hz 107 kPa 
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Figure 23 

Series 4b Increasing Frequency 
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HELIUM EXPERIMENTAL 1077 Hz 107 kPa 
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Figure 24 

Series 4c Increasing Frequency 
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Figure 25 

Series 4d Increasing Frequency 
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HELIUM EXPERIMENTAL 1860 Hz 167 kPa 
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Figure 26 

Series 4e Increasing Frequency 
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(KX in multiples of nTT) and degradation at the pressure nodes 
(KX at ( 2n+l ) ( 7T/2 ) ) , Similar skewing to what he showed was also 
evident • 

A second feature present in all of the series was an 
asymmetry in the magnitude of the maximum temperature 
differences. The magnitude of the positive temperature 
difference (Tj, ~ ) was always less than that of the negative 

(T(- - T„ ) difference. 

The shape of the curves also exhibited interesting 
properties. The curve shape was remarkably similar from cycle to 
cycle on the same side of the At axis, but varied greatly from 
above to below the axis. 

Examination of the maximum temperature differences vs 
pressure amplitude (Figure 27) reveals the rate at which the 
temperature difference falls off as the acoustic and mean 
pressure increases. Several of the data sets start to show a 
saturation of the effect at high acoustic and mean or static 
pressures . 
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Figure 27 

Maximum Temperature Difference vs Pressure Amplitude 
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V. SUMMARY/RECOMMENDATIONS 



The main accomplishment of this thesis was to put together a 
system enabling computerized measurement of temperature gradients 
developed across short plates by thermoacoustic heat transport. 
Measurements were made of the gradients under different static 
and acoustic pressures, different frequencies, and different 
gases. Comparison with theoretical values showed results similar 
to those found by Muzzerall and further showed a gradual 
degradation of agreement from low to high acoustic pressures. 

The thermoacoustic effect becomes saturated at the higher values 
of acoustic pressure amplitude. Also, theoretical values and 
measured data diverge most severely in the area of pressure 
nodes. The causes of both of these effects bear further 
investigation. 
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